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G ' Using the fermion-spin transformation to implement spin-charge separation of 



constrained electrons, a model of two t — J chains with interchain single-electron 

hopping is studied by abelian bosonization. After spin-charge decoupling the charge 
(N 

^ ' dynamics can be trivially solved, while the spin dynamics is determined by a strong- 

'nT ■ 

(<— ^ ' coupling fixed point where the correlation functions can be calculated explicitly. 

^^ ' This is a generalization of the Luther-Emery line for two-coupled t — J chains. The 

■^ 



interchain single-electron hopping changes the asymptotic behavior of the interchain 
spin-spin correlation functions and the electron Green function, but their exponents 



'^ ■ are independent of the coupling strength. 

o 



. c^ ; PACS numbers: 75.10.Lp, 75.10.Jm 



An important issue of current interest is whether the pecuhar properties of one- 
dimensional (ID) Luttinger hquid (LL) [|T|, [|2[ will survive in two- and three-dimensions. The 
renormalization group (RG) studies seemed to indicate an instability of LL behavior with 
respect to the interchain single-electron hopping (SEH) t± 0. However, Anderson suggested 
that one should treat the intrachain correlations exactly including the spin- charge separation 
before switching on SEH. Using the asymptotic Green functions in ID for a finite Hubbard 
t/, he argued that SEH is an irrelevant variable and named this property as "confinement" 



of the ID Hubbard model 0. His idea has stimulated several further studies [g- H, most 
of which did not confirm his conjecture in the strict sense. SEH is indeed renormalized 
to zero, but e-e or e-h pair hopping is generated, which drives the coupled chains towards 
a strong-coupling fixed point corresponding to superconducting or density-wave states. In 
fact, this type of instability was studied earlier in connection with organic superconductors 



TT| . Nevertheless, this result is not convincing because the validity of the perturbative RG 
at strong-coupling fixed point with large U is questionable, as in the single-impurity Kondo 
problem. The Kondo physics is determined by the Wilson strong-coupling fixed point []T2 . 



The poor- man's scaling []T^ correctly directs the RG flow towards it, but the calculation can 



not be justified by itself |T^ . 

In this paper, we consider two coupled t — J chains, using a fermion-spin transformation, 
proposed recently by Feng et al. ||15[|, where the charge degrees of freedom are described by 
spinless fermions, while the spin degrees of freedom are represented by hard-core bosons, 
which in turn, can be expressed as another type of spinless fermions via Jordan- Wigner 
transformation. The on-site local constraint for single occupancy is satisfied even in the 
mean-field approximation (MFA) and the sum rule for physical electrons is obeyed. We 



combine this transformation with the abelian bosonization technique |jT|], |[T6| to consider the 
effect of SEH on the correlation functions. After spin-charge decoupling the charge dynamics 
can be solved trivially, while the spin dynamics can be mapped into noninteracting spinless 
fermions. This strong-coupling fixed point is similar to the Luther-Emery line of the single 
chain problem with back scattering fl^ . We confirm that the spin-charge separation by itself 



does not produce Anderson confinement |^, 0. Moreover, SEH changes the asymptotic 
behavior of the interchain spin-spin correlation functions and the electron Green function, 
but their exponents are independent of the coupling strength t±. 
We consider two coupled t — J chains 



i,(7 i,(T 



with local constraint J2a C'i>C'i,o-<l- Here C*} j ,^(C2^j ,^) creates an electron with spin a at site 
i on chain 1 (2), and Si^i{S2^i) is the corresponding electron spin operator; ty is the intrachain 
hopping and /i is the chemical potential. The fermion-spin transformation of constrained 



electrons [|15| 



can implement the spin-charge separation without additional constraints. Here Oj and a\ 
are "ho Ion" (or "electron" in the particle representation) operators, represented by spinless 
fermions. Sf and Sf are spinous or pseudo-spin operators represented by CP^ hard-core 
bosons, different from the electron spin operators in Eq.(l). P is a projection operator 
removing the extra degrees of freedom in the CP^ representation. The anticommutation 
relations for constrained fermions Ci^^ are strictly preserved. Moreover, the local constraint 
is satisfied exactly. However, the projection operator P is cumbersome to handle and in 
many cases, for example, MFA, we can drop it, with very good results JT3[ . 
To establish notations, consider first a single t — J chain 

i i 

-F2J^(aJai)S'iS'i+i(a|+iaj+i). (2) 

i 

As shown in |T^, using the above fermion-spin representation, the Jordan- Wigner transfor- 
mation St = //e^"S'<> -^'^^S ^r = {StV, Sf = fjfi - i and MFA, one finds the ground state 



energy and gapless spinon and holon spectra, in good agreement with the exact solution [|18 . 



However, to obtain correct exponents for correlation functions, one has to go beyond the 
MFA, taking into account holon-spinon interactions. Following Weng et al. [|I^, this can be 
done by "squeezing out" holes from the spin chain, i.e., to replace ajaj+i(/j /j+i + /i/j+i) 
by ajai^i wherever there is a hole at site i and introducing the "string operators" which in 
our case are given by 



an 



In the resulting Hamiltonian, the "holon" part is free and can be easily bosonized, while the 
spinon part is an antiferromagnetic Heisenberg spin-1/2 chain, which can also be bosonized 



and reduced to a standard 1+1 quantum sine-Gordon (SG) model ^0\ 



Hu 



/..(-^n^ + ^(V.)^-^co.v/l6^,), (3) 



where a is an ultroviolate cut-off, while the boson field ip describes the low-energy excita- 
tions of spinous, n is its conjugate momentum with a commutation relation [(f{x), n(a;')] = 
i6{x - x'). The spinon velocity is Vg = 2J (1 - 6^ - (^^^) Jl + ^, with 6 as doping 
concentration. The parameter determining the exponent of the spinon correlation function 
is -ft's = (1 + -)~^/^, which should be independent of 6, and our result for Kg is slightly away 
from the exact value derived for half- filling pO|, [^. In principle, the abelian bosoniza- 



tion is exact only at J^/Jl^O for the Heisenberg spin-1/2 chain [|T6|. However, the exact 
Bethe-ansatz solution does not show any singularities for —1 < Jz/J± < 1, and the isotropic 



antiferromagnetic coupling Jz/J± = 1 is described by this fixed point [^. On the other 
hand, the 1+1 SG model with AcosjSip has only one weak-coupling fixed point for small A 
at /5^ = Stt [E^. Thus we can associate the cosine interaction (Rl) with this fixed point of 



the SG model in order to rectify Kg to be 1/2 after rescaling H -^ v'^H, ip —>■ -^. Since 
the fixed point of SG Hamiltonian under RG for (3 > 8tt corresponds to the vanishing of the 
cosine term [^, we can easily calculate the asymptotic behavior of the spin-spin correlation 



functions and the electron Green functions, in good agreement with exact results [T^, e.g. 



S{x, - x„ t) cos2kF{x, X,) _^ 



(TQ,.(t)C],,(0)) 



e 



(5) 



[{Xi - XjY - {Vhty] 16 [{Xi - Xj) - (Vst)] 2 [(Xj - Xj) - (Vht)] 

where the holon velocity is the exact value Vh = 2t^\sin6TT, and kp = |(1 — 6). 

Now consider two coupled chains and use the MFA to decouple the interchain holon- 
spinon interaction. The Hamiltonian (1) is reduced to the following form H = H^ + Hs, 
and 

Hh = -t\\ ^(al^jai,i+i + 4^ja2,i+i + h.c.) - t±r]i Xl(4,ia2,i + h.c.) 

i i 

-/i ^(4 jOi^i + alia2,i), (6) 

i 

Hs = 2r^f Y.iSlA,^+l + S2A,^+l) - t±V2 T.iStiS2, + h.c), (7) 

i i 

where we have defined two MF order parameters rji and ri2. 

The holon Hamiltonian is trivially diagonalized by introducing A^ = 4s(fli,fc + 02,fc) and 
Bk = —-j^{ai^k — o.2,k) with excitation energies e^ = —2t\\cosk — t± and ef = —2t\\cosk + t±, 
respectively, where we assume rji ^ 1, as will be confirmed later. The SEH splits the 
original holon excitation spectrum by 2t± and in the low doping case for a finite value 
t± > t\\{^~ cos26it), only the upper band has vacancies and the lower band is fully occupied. 
The above condition on t± is usually satisfied. Thus, it is easy to find the self-consistent 
value ri2 = —{1 — 6), a.s well as the interchain holon correlation functions using the abelian 
bosonization technique |]I| 



/' 



{x, - x)r - {v,tY -, (8) 



(e-lE,<.4,(*K,^Wa,^^(t)4.(0)e^5E,<,«l(0KK0)^ 



e 



tkp\Xi X .) 



{xi - x\y - {vhty "" {xi - x'-) - (vht) 



(9) 



where kf = kp — ^- Due to the single occupancy constraint, the situation here is simpler 
than the weak-coupling case where both bands have to be considered |^- [ p!0[] . 
The spinon part can be reduced to the following form: 

+^V^cos(V^(s5,-s5,))], (10) 

where (p{x) is the dual field of (p{x), and is defined by || = 11 and H = — |^. Note the 
difference of (0) from (|^), where no dual fields are involved. Introducing symmetric and 
antisymmetric combinations of spinon fields: ips = tToIv^i + </'2), fA = "^(v'l ~f2), (|10D can 
be rewritten as: Hg = H^ + Hf, where 

Hf^><^n| + ^(v..)^), (11) 



'VsKs 2 , ^s / ^2 , (1 -^)^-L 



Hf= dxl 4-^ni + -iriV^AY + ' , ^^ cos{V27TipA)]. (12) 



2 "" 2K,^'^ ' (na 
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The symmetric part H^ is a LL with the same parameters Vg and Kg as for a single chain. 
As for the antisymmetric part, a rescaling and use of self-duality of the non- interacting part 
lead to the following 1+1 quantum SG model 

'vgr-P^o , ^ ,0, (1 — (5)i± / /27r 
^mU (V^.)'] + ^^ca,( J- 



ff." = / -iM ^|n^ + (V0.)'] + +rr|^"»(i/l^45^) I ■ (13) 



With the corrected -ft's = |, the coupling strength of the SG model is (? = Atx. This is 
nothing but a free massive Thirring model [^] with a mass gap in the excitation spectrum 



Ag K, 2(1 — (5)t_L. It is known that t_L is a relevant variable in the range < /?^ < Svr. Some 



years ago, Haldane [23| considered the renormalization of the Bethe ansatz equations for 



the massive Thirring model, equivalent to the 1+1 SG model [Q. He extracted a quantum 
fiuctuation parameter that controlled the correlation functions of this model, and found 
that at (3'^ = Att the renormalization of the model stops and it corresponds to a free spinless 
fermion field. This means that (3"^ = Air is just the strong-coupling fixed point , analogous 
to the Toulouse limit of the single-impurity Kondo problem ^^ . It is the fixed point that 



6 



controls the properties of this model in the whole region < (3"^ < Stt. It is remarkable that 
after correcting the Ks value using the Bethe ansatz solution for the single chain, we end 
up exactly at this fixed point for two coupled t-J chains. If we did not rectify the parameter 
Kg in the absence of the SEH, the coupling strength of the above SG model would be at 
(/5 )^ < Air, and it should be renormalized to strong-coupling fixed point /5^ = 4it, while the 
parameter Kg is renormalized to |. In the end, the spinon correlation functions of the single 
t — J chain could still have correct asymptotic behavior. 

In the weak coupling approach, it was also found that the ground state of a single 
chain is unstable with respect to SEH and the two-coupled chains are driven to a strong- 
coupling fixed point corresponding to /5^ = Air with opening a gap in one of the spinon 
excitation spectra |Q, which is equivalent to the Luther- Emery line for a single chain with 
back scattering. However, for a finite t±, the spin and charge degrees of freedom are still 
coupled, and the renormalizaion process cannot be carried out in the perturbative approach. 
Therefore, it is not possible to calculate the correlation functions at this strong-coupling 
fixed point . On the contrary, within our strong-coupling approach, we end up exactly 
at this strong-coupling fixed point and the spinon excitations reduce to two modes, one 
is the original LL branch for the single chain, while the other is a massive free fermion 
branch, corresponding to the soliton gas of the quantum SG model. After convolution, the 
interchain spin-spin correlation function and electron Green function for two-coupled chains 
can be calculated as 

S{x, - x' t) cos2kF{x, - x/) _ 

[(x, - x/y - Kt)2]T6 [{x, - x-Y - {vgtYY 



{TC,,At)cLM) 



1 

e""F 



fp \'^i "^ j ) 



(Xi - x'jY - {Vhtf "^ {Xi - x'j) - (Vst) * [{Xi - x/) - (Vht)] 2 



(15) 



The parameter r]i ^ 1, a.s mentioned earlier. As compared with (^) and (|^) for a single 
t — J chain, the SEH has generated new exponents, independent of t±. The singularity 
is weaker due to the presence of a gap in one of the excitation branches. The "spinon" 



exponent is — | instead of —1 for the spin-spin correlation function, while it is — ;| instead of 
— I for the electron Green function. The singularity due to holons is also weakened because 
of the hybridization of two chains and the single occupancy constraint. The gap in the 
excitation spectrum for the antisymmetric spinon field (f) leads to an exponential decay for 
(0(a;, t)(/)^(0, 0)), but a constant contribution to the correlation functions for the fermionic 
fields. 

To summarize, we have found strong-coupling fixed point controlling the behavior of 
two coupled t — J chains and have calculated explicitly the interchain spin-spin correlation 
function and the electron Green function. Our work has reconfirmed some results of the 
previous weak-coupling studies, namely, (i) the spin-charge separation does not produce by 
itself Anderson confinement; (ii) the exponents of the interchain correlation functions are 
changed due to the presence SEH P, 0. However, there are significant differences between 
these two approaches, (i) The weak-coupling approach indicates the existence of a strong- 
coupling fixed point , but can not provide a valid calculation scheme at that fixed point, 
(ii) The spin-charge separation for a single chain in the weak-coupling sense in general does 
not guarantee the spin-charge separation for two coupled chains (except for a special ge- 
ology model [Q), so it is not possible to calculate explicitly the correlation functions. The 
situation here is similar to the single chain problem. The correlation exponents depend on 
the interaction strength in the weak-coupling limit, while it is independent of interaction 
strength in the strong-coupling limit; the spin-charge separation is valid in the sense of 
" almost complete factorization" of the wave function in the large-U limit [l^ . Here we use 
the spin-charge separation in the same sense and it can be thus justified. 
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